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CHAPTER  1 


INTRODUCTION 


The  object  of  this  study  Is  to  develop  an  efficient  numerical 
solution  for  the  high-frequency  radiation  patterns  of  a  spherlod- 
mounted  antenna.  The  Geometrical  Theory  of  Diffraction  (GTD)  Is 
employed  for  the  pattern  calculation  of  which  the  creeping  wave  [1] 
solution  In  the  shadow  region  Is  of  particular  Interest  here. 

According  to  the  generalized  Fermat's  principle,  a  ray  emanating  from  a 
source,  which  Is  located  on  the  surface,  follows  a  geodesic  path  on  the 
surface  and  continually  sheds  energy  Into  the  shadow  region.  Such  a 
creeping  wave  mechanism  Is  Illustrated  In  Figure  1(a),  from  which  It 
can  be  seen  that  a  ray  traverses  from  the  source  point  Q’  to  the 
diffraction  point  Q,  and  then  propagates  along  the  geodesic  tangent  at 
Q  toward  the  observation  point  P?.  Since  most  parameters  Involved  In 
the  pattern  calculation  are  related  to  the  ray  path,  the  geodesic 
solution  on  the  surface  constitutes  a  major  portion  of  this  study. 


LIT  REGION 


Figure  1.  Ray  paths  In  the  shadow  and  lit  regions. 


The  spheroid  is  studied  because  of  Its  three  dimensional  nature 
in  that  It  can  be  used  to  simulate  an  aircraft  or  missile  fuselage. 

For  a  flush-mounted  airborne  antenna,  the  fuselage  has  a  dominant 
effect  on  the  resulting  radiation  pattern.  Previously  [2]  the  fuselage 
was  simulated  by  an  elliptic  cylinder;  however,  the  cylinder  model 
could  not  predict  the  pattern  close  to  the  nose  or  tall  sectors,  where 
the  deviation  from  the  physical  situation  becomes  very  prominent. 
Previous  attempts  [2],  using  the  spheroid  model,  employed  calculus  of 
variations  or  tensor  analysis  to  calculate  the  geodesic  paths.  The 
resulting  algorithms  were  very  time-consuming  and  virtually  impractical 
to  use.  Therefore,  a  perturbation  method  Is  presented  here  In  an 
attempt  to  combine  the  simplicity  of  the  cylinder  or  cone  solution  with 
the  generality  of  the  prolate  spheroid. 

Before  examining  our  perturbation  approach.  It  Is  appropriate  to 
study  the  basic  radiation  mechanisms  In  some  detail  In  that  they 
suggest  the  present  solution.  The  ray  analysis  in  the  lit  region  can  be 
rather  straightforwardly  analyzed  for  ^ny  convex  surface;  however,  the 
geodesics  which  are  necessary  to  treat  the  shadow  region  dictate  what 
surfaces  can  be  efficiently  analyzed.  In  the  shadow  region,  the  energy 
propagates  outward  from  the  source  along  the  geodesic  paths.  As  the 
energy  flows  around  the  surface.  It  Is  continuously  diffracted  along 
the  geodesic  tangent  toward  the  field  point  such  that  the  significant 
effect  of  the  surface  Is  associated  with  a  region  around  the  source. 

In  fact,  for  a  prolate  spheroid,  the  significant  portion  of  the  surface, 
which  Is  associated  with  the  dominant  energy,  may  look  as  shown  in 
Figure  2.  Note  that  this  region  can  be  specified  by  following  the 
various  geodesic  paths  until  the  radiation  level  along  a  given  path 
becomes  insignificant,  i.e.,  more  than  40  db  below  the  pattern 
maximum.  With  this  In  mind.  It  Is  clear  that  one  could  represent  the 
prolate  spheroid  by  a  structure  which  simulates  the  circular  cross- 
section  completely;  however,  the  profile  (I.e.,  the  elliptic  curve) 
could  be  approximated  by  a  simpler  shape  In  that  the  significant  energy 
region  does  not  cover  a  large  portion  of  the  profile  shape. 

Consequently,  a  right  circular  cone  perturbation  Is  chosen  here  to 
simulate  the  spheroid  which  In  turn  can  be  used  to  model  a  fuselage. 

This  perturbation  model  Is  Illustrated  In  Figure  3  for  a  source 
located  near  one  end  of  the  spheroid.  Note  that  If  the  antenna  is 
placed  at  the  center  of  the  spheroid,  the  cone  actually  becomes  a 
right  circular  cylinder. 

In  that  the  cone  Is  a  developed  surface,  one  can  unfold  the  cone 
such  that  a  planar  structure  results.  The  geodesics  associated  with 
the  cone  are,  then,  straight  lines  on  this  planar  structure.  In  order 
to  allow  for  a  geodesic  solution  between  the  simplicity  of  the  cone 
and  the  rigor  of  the  spheroid,  one  can  perturb  the  cone  by  bending  It 
along  Its  generator  as  Illustrated  In  Figure  3(b).  In  that  a 
perturbation  technique  is  employed,  the  geodesic  paths  for  the  cone  are 
simply  modified  such  that  the  solution  for  the  spheroid  is  basically 
straightforward  and  Is  the  basis  for  this  thesis,  It  Is  obvious  that 
one  cannot  use  this  perturbation  If  significant  energy  propagates  great 
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distances  along  the  surface.  However,  as  discussed  above,  the  energy 
which  propagates  great  distances  along  the  spheroid  generator  becomes 
Insignificant  In  magnitude  such  that  one  need  not  solve  for  the  true 
geodesic  paths  outside  the  region  shown  In  Figure  3.  The  simplicity 
of  these  perturbed  geodesic  paths  allows  one  to  very  efficiently 
determine  the  significant  ray  paths  on  the  prolate  spheroid. 

It  Is  noted  that,  since  GTD  Is  a  high-frequency  method,  the  lower 
frequency  limit  of  this  solution  Is  dictated  by  the  electrical 
dimensions  of  the  spheroid,  l.e.,  the  semi-minor  axis  of  the  spheroid 
Is  required  to  be  at  least  a  wavelength. 

The  curved  surface  diffraction  solutions  employed  here  are 
briefly  discussed  In  Chapter  II.  The  parameters  needed  for  this  GTD 
solution  are  described  In  Chapter  III  In  terms  of  the  cylinder  and  eoi 
perturbations  used  to  simulate  the  geodesic  paths  on  the  spheroid. 
Analytic  results  are  compared  with  various  measurements  In  Chapter  IV 
Finally,  a  summary  of  this  study  Is  made  in  Chapter  V. 


CHAPTER  II 


THEORETICAL  BACKGROUND 


A.  INTRODUCTION 

The  radiated  field  of  a  spheroid-mounted  antenna  is  analysed 
using  the  Geometrical  Theory  of  Diffraction  (GTD).  The  surface  is 
assumed  to  be  perfectly  conducting,  and  the  surrounding  medium  is  free 
space.  An  exp  (jwt)  time  dependence  is  understood  and  suppressed  in  the 
following  formulations. 

Consider  an  infinitesimal,  magnetic  current  moment  dP’  (Q')  or  an 
electric  current  moment  dFe(Q')  located  on  a  perfectly  conducting 
convex  surface  as  shown  In  Figure  1;  the  sources 

dFm(Q')  -  F(Q' )xn'da' 


and  dFe(Q')  =  I(*  '  )d«,'  n' 

pertain  to  the  aperture  and  monopole  type  excitations  with 
E(Q')  =  electric  field  at  Q' , 
n'  =  outward  unit  surface  normal  at  Q' , 
da'  =  area  element  at  Q' , 

I(fc')  =  electric  current  distribution  on  the  monopole,  and 

*'  =  distance  parameter  along  the  monopole. 

According  to  geometrical  optics,  the  space  surrounding  the  source  is 
divided  Into  an  illuminated  and  shadow  region  by  a  plane  tangent  to  the 
surface  at  Q' .  This  plane  is  referred  to  as  a  shadow  boundary.  Since 
the  field  in  the  deep  lit  region  is  essentially  that  obtained  from 
geometrical  optics,  and  the  field  In  the  deep  shadow  region  is 
relatively  weak,  the  solution  for  the  transition  region  adjacent  to  the 
shadow  boundary  is  of  more  Interest  and  discussed  below. 


B.  SHADOW  REGION 

The  creeping  wave  mechanism  in  the  shadow  region  is  Illustrated 
in  Figure  4.  From  the  generalized  Fermat's  principle,  a  ray  emanating 


from  the  source  dPrg(Q')  at  Q'  traverses  a  geodesic  oath  TO  on  the 

surface,  and  propagate^  along  the  geodesic  -?nqent  at  Q  toward  the  field 
point  P  .  The  field  dEg  at  P$  can  be  expressed  in  terms  of  the  field  at 

a  reference  point  PQ  by  [3] 


r  pdpd~  ' 

dljs(PsMFs(P0)  4 

v  (pi+soHp2  so} 


■2 


where 

and  p«  are  the  principal  radii  of  curvatures  of  the  wavefront 
1  2  ^  P0; 

and  0[m‘2,m"3,. . .1  are  the  higher  order  terms. 


From  Fiqure  4,  it  is  seen  that  if  the  reference  point  P0  is  moved  to  the 

3  d  d 

curved  surface  diffraction  point  Q,  then  p-j->0,  P2'*pc*  an<*  s0‘*'s-  ^nce 

dEm(Ps)  is  independent  of  the  reference  point  PQ,  It  follows  that 

11m  /p?  dFm(P0)  =  r-(Q',Q)  ;  (2) 

P0-*Q’  4 

then, 

d^(PS)^(OMl)  jj^«7e'3kS  '  131 

Furthermore,  T^Q'.Q)  can  be  related  to  the  source  strength  dP^  at  Q' 
by  e 

^(Q’,0)  =  dPg  (Q')-f^(Q',Q)  (4> 

where  f|(Q',Q)  is  given  by  [3] 
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Tm(Q',Q)  «  ^|^£b,nT1  (Q * ) H+t ' bT2(Q ' )S+6 » 6t3CQ '  JS+t'n^CQ'  )H] 


e-jkt  <**o 

Pg(Q') 

fe(Q'.Q)  »  ^2-  [n * nT 5 ( Q ' ) H+n 1  bT 6 ( Q ' )S]  e"Jkt 

l~\~r pa(Q)  J/6 


(5) 


(6) 


*  A  A  AAA 

Here  (t‘ ,n',b')  and  (t,n,b)  are  the  tangent,  normal  and  bl normal  unit 
vectors  to  the  surface  at  the  source  point  (Q* J^and^ diffraction  point 
(Q),  respectively.  As  seen  from  Figure  4,  t  x  n  e  b  and  t'  x  n'  *  6*. 
The  quantities  T« (Q‘ ) »• • • »Tg(Q' )  are  the  torsion  factors  at  Q'  and  are 
given  in  Table  1 !  Also, 


h  *  g(0 

(7) 

**mM<*> 

(8) 

with 

g(e) 


dr 

°°exp(- j2ir/3) 


W^(t) 


and 


(9) 


g(0 


dTex£tMl 

«oexp(-j2Tr/3)  w2u; 


(10) 


which  are  known  as  the  acoustic  hard  and  soft  Fock  functions.  The  Fock 
type  Airy  function  Is  given  by 


*2^  a  7?  J  dt-exp  (xt-t3/3) 

2  ^  -exp(J2w/3) 
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(11 


that  appear  1 
Geometry  [4]. 


and  w£(t)  Is  the  derivative  of  w2(t)  with  respect  to  r.  The  Fock  para¬ 
meter  g  for  the  shadow  region  Is  given  by  [3] 


02) 


03) 


Here  oq(t')  Is  the  surface  radius  of  curvature  along  the  ray  path  at  t'. 
The  width  of  the  surface  ray  tube  at  Q,  dn(Q)t  Is  given  by 


dn(Q) 


pc<4 


(14) 


The  parameters  Z  and  t  are  defined  as  the  free  space  wave  Impedance 
and  geodesic  arc  legth  from  Q*  to  Q,  respectively. 

Combining  Equations  (3)-(14),  the  rt  and  b  directed  components  of 
dt|(Ps)  ere  given  by  [3] 


(a)  dPm((|*)  case: 


n  Ik  -  4k*rP0(Q,)/V6pi:r 

<<M  *  ^^'b' jar  Jan^sT 

e‘Jks+0[m“2] 


05) 


d0Ps)  *  C(A-b)T0Sf(dfc.t*)S]  e‘jktC^l]“1/6 


dip 


__0  /  1  .-JkSxnrm'2 

dip  Js(pVs)  e  ^Cm  ,m  ] 


(16) 
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(b)  dPe(Q')  case: 


<<Ps>  - 

e‘-iks+0[m“?l 


0$T]  l~®  sTp^+sT 


(17) 


b  -Jkz 
d£e<ps)  *  ”4? 


o  -ikt  pa(Q'>  -V6  fK 

^  dPe(Q')T0Se  '“liV-l  ,/aj2 


3?”  /s(p>s) 


e"iks+C[m'2] 


08) 


where  To=T(Q')p  (Q')  with  T(Q')  being  the  surface  torsion  at  the  source 
location  (referyto  Table  1). 


C.  LIT  REGION 


which 

point 

point 


From  geometrical  optics,  the  source  dRn(Q')  at  Q'  excites  waves 
propagate  along  straight  line  ray  paths  from  the  source  to  field 
In  the  lit  region.  As  shown  In  Figure  5,  the  field  dEro(P, )  at 
PL  Is  expressed  by  e  L 


dFg(PLhdtjg(P0) 


(19) 


Since  Q'  Is  the  only  caustic  of  the  Incident  rays,  the  principal  radii 
i  i 

of  curvature  p,  and  p9  associated  with  the  Incident  wavefront  at  P„  are 
1  4  2.  ,  _  o 

Identical,  1.  e.,  pj*©^*  p-j.  Furthermore,  dEm(PL)  is  independent  of 
the  reference  point  P  .  If  PQ  Is  chosen  to  be  at  Q',  it  follows  that 


11m  pi 

V*' 

p1"*0 

Vs 


dEg(Po) 


i 

Ig 


(20) 


12 


I  1 

I  i: 


i: 


i: 
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A  A 

N£^i; 


»  •  S'  cotff1  * 

bj ■  *4  *f*b  ■  tun 


(b) 


Figure  5-  Ray  tube  and  ray  coordinates. 


should  exist.  Thus,  can  be  related  to  dPjg(Q')  by  [3] 


-it 

i-m 


W& 


Equations  (19)-(21)  are,  then,  contained  to  yield 


d^PLl^^g' )  ►0[m^2,m"13, . . .] 


iL  A  A  A  A  A  A  A  A 

Tm  '  iT  Cb;nA+t;bB+bJbC+tJnD] 
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(21) 


(22) 

(23) 


-*,£H  i  '  ^  !£"T  e*'  ft  t T 

1>miivS>SSB83bMK 


(24) 


e  4n 


[n'nM+n'bN] 


with  A,  B,  C,  D,  M,  and  N  are  defined  In  Table  2.  Note  that  df^(PL)  Is 
decoupled  Into  n  and  b  components  as  follows: 


(a)  dPm(Q')  case: 


dEm(pL)  =  =£  C(dPm* b ' ) (HJt+To2Fc°s01 )  +  (dFm-t'  ^cose1  ] 


.-jks  2 

+  0[mt2] 


(25) 


<(pl)  * 


“jkS  o  o 

— T“+  oCnV  A  1 


(26) 


(b)  dF.(Q')  case: 


—  Ak7 

dEe(PL)  *  “THT  dpe(Q')s1neTtH£+Tc^cos0l]  +  °K^ 


1  rut . t  2c____1t  e 


jks 


-2- 


(27) 


-jk7 


jks 


dEe(PL>  *  -T71  4pe(Q')s'"«'V  H-  * 


(28) 


where 


»  -j^/3 
w  *  g(cA)e 

S,’S^T^t)e 


(29) 

(30) 
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I 


ei  *  -^(Q'Jcos0 


mf<Q‘)  “ - g(qi-4  m 

1  (l+VcosV)1'3 


The  angle  Is  defined  by  n » . £  ■  cosfl*  as  shewn  In  Flqure  5, 


r  _  Sl-H*cose 

- 5 — TT 

1+T  * cos 


as  defined  In  Table  II. 
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Also 


D.  PATTERN  FACTORS 

The  solutions  for  short  magnetic  or  electric  dipoles  have  been 
given  In  part  (B)  and  (C).  One  approach  to  analyze  an  extended  aperture 
or  linear  antenna  problem  Is  to  Integrate  the  above  solutions  over  the 
source  distribution.  If  It  Is  known.  This  Is  an  application  of  the 
superposition  theorem,  and  one  approximates  the  source  distribution  by 
an  array  of  short  magnetic  (or  electric)  dipoles  on  the  conducting 
surface.  This  Is  an  accurate  solution,  however,  rather  tedious.  A 

more  efficient  approach  Is  to  modify  dP^(Q')  as  shown  In  Reference  5 
such  that 


(a)  in  the  shadow  region: 


-  p  26 

Km  ~  ~ 


A  A 


cos(j£(PB.t')) 


'-irP? 


s1n(!£  -b* ) 


III 


p  .5  • 

_  2  r>' 


*  n' [l-cos(kl)l 


(35) 


(b)  in  the  lit  region: 


F  *  P  — 

r  m  *  rr» 


cos  (?p«  1  ne 1  ( P(1*  t ' ) )  j 

;  s1n(|^s1ne1(Pn.*b')) 

Mj-sine^Vt'))2! 

-^slneV-'b’) 

(36) 


1?  -  cos(kLn' »s)-cos(kL) 

r  “  M  *  *  O 

e  l-(n' *s)Z 


(37) 


Here  P„  *  unit  vector  In  the  direction  of  magnetic  current  moment. 

Ml 

A,B  *»  the  length  of  the  short  and  long  sides  of  the  slot,  and 

L  =  the  length  of  the  monopole. 

It  Is  noted  that  L  Is  not  to  exceed  a  quarter  wavelength  for  the 
solutions  to  be  valid. 
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E.  SPHEROID  SURFACE  PARAMETERS 

The  formulations  In  sections  (B)-(D)  are  used  to  solve  for  the 
radiated  fields  of  antennas  mounted  on  a  prolate  spheroid.  Using  the 
spheroid  geometry  shown  In  Figure  6,  the  surface  Is  defined  by 


R(e,$)  *  R(e)s1necos$x+R(e)s1nes1n$y  +R(e)cosez 

with 

R(e)  3  - - -ac - 

/a^cos2e+c2s1n2e 


or, 

“*■  A  A 

R(v,$)  *  a  cosvcos$x  +  a  cos\>s1n<j>y  +  c  slnvz 

wl  th 


tanv  a 


c  cose 
a'  sine 


(38) 


(39) 


(40) 


(41) 


The  "v"  parameter  Is  Introduced  because  of  its  convenience  in 
analyzing  elliptic  geometries. 

Considering  a  ray  which  propagates  along  a  geodesic  path  Q'Q  on 
the.,  spheroidal  surface  as  shown  In  Figure  7,  the  three  unit  vectors 
t,  n  and  6  are,  as  defined  earlier,  the  geodesic  tangent,  outward 
surface  normal  and  blnormal  at  any  point  along  the  geodesic  path.  The 
outward  surface  unit  normal  (n)  Is  obtained  from 


“  _  R&xRv 

j  $$xftv  | 


(42) 


where 


|5-Sv- 

3v 


-a  slnvcos^x  -  a  slnvsln^y  +  c  cosvz 


and 


8-* 


-a  cosvslntx  +  a  cosvcos$.y 
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Geometry  of  a  spheroid. 


Figure  7.  Geodesic  path  on  a  spheroid. 
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Then, 


n 


A  A  A 

_  c  cosycositix  +  c  cosvslmy  +  a  sinvz 


I 


2  2  2  2 
a  sin  v  +  c  cos  v 


(43) 


The  normal  curvatures  *  on  the  surface  are  evaluated  by  introducing 
the  first  and  second  fundamental  forms  of  differential  geometry  [4] 
such  that 

_  L<*t>2  +  2Md<t>dv  +  Ndv2 

k - - - -  (44 

n  Edr  +  2Fd$dv  +  Gdv^ 


where 


L  ■  V"1  M  =  V*'  N  =  V 


E  =  W*  F  =  V*  .  6  =  t  -t 
$  <p  v  v 


N 

v 


and 


It  can  be  shown  that 

=  -a  cosvcostx  -  a  cosvsin<t>y 

->•  A  A 

R^v  =  a  slnvsin^x  -  a  sinvcos^y  ,  and 

"ft  =  -a  cos\>cos<f>x  -  a  cosvsirtfy  -  c  sinvz 
w 

After  some  algebraic  manipulation,  one  obtains 

2  r 2  2  2  2 

L  =  ac  cos  v//a  sin  v  +  c  cos  v 


M  =  0 
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N  *  ac//a2s1n2v  +  c2cos2v 


2  2 
E  =  a  cos  v 


F  =  0 


G  =  a2s1n2v  +  c2cos2v 


The  two  principal  curvatures  on  the  surface  are,  then,  given  by 


*  L/E  = 


/  2  2  2  2 
a  sin  v  +  c  cos  \j 


(45) 


and 


N/G  = 


ac 


(a2s1n2v  +  c2cos2\>)^ 


(46) 


which  can,  also,  be  expressed  as 


_  c  fa2cos2e  +  c2s1n2e 
1  a  / a4cos2e  +  c4s1n2e 


(47) 


and 


_  ac(a2cos2e  +  c2s1n28)^2 
(a*cos2e  +  cSln2er^ 


(48) 


It  Is  noticed  that  R,  =  —  and  R,  *  —  as  found  In  Table  I.  Since 
F=M=*0,  the  directions  R.  and  ft  are  In  the  principal  surface 

9  V  A  A 

directions  [4]  which  are  expressed  by  TpT2  and  are  associated  with 
and  k?»  respectively,  as  shown  in  Figure  7.  If  0  denotes  the  angle 

“  A  A 

between  t  and  ,  then 


A  A 


t  *  r-jcose  +  T2Slng 


(49) 
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From  Euler's  theorem,  the  normal  curvature  along  the  geodesic  path  Is 
specified  by 


2  2 
?g  a  Kjcos  6  +  <2sin  e 


(50) 


with  the  radius  of  curvature  (pg)  being  l/r^. 

The  torsion  term  (T  )  Introduced  In  Section  (B)  Is  given  by 


T0  ■  T-Pg 

where  the  surface  torsion  Is  given  by 

t  -  s1n2$  /  ^ 

T - (k2-k}) 

with  < i  and  k.2  being  defined  In  Equations  (45)-(48). 


(51) 


(52) 


CHAPTER  III 

NUMERICAL  TECHNIQUE  AND  PERTURBATION  METHOD 


A.  INTRODUCTION 

It  Is  seen  that,  for  an  antenna  mounted  on  a  spheroid,  the 
solution  In  the  lit  region  Is  relatively  straightforward.  On  the  other 
hand,  the  geodesic  paths  associated  with  the  GTD  solution  In  the  shadow 
region  are  extremely  complex.  Previous  studies  [2]  employed  calculus 
of  variations  or  tensor  analysis  to  calculate  the  oeodeslc  Daths  which 
resulted  In  very  Inefficient  solutions.  An  efficient  numerical 
approach  Is  examined  In  this  chapter  with  the  spheroid  simulated  by  a 
perturbed  cone  or  cylinder  model.  In  that  the  cone  and  cylinder  are 
developed  surfaces,  geodesics  are,  then  efficient  to  solve.  Given  a 
radiation  direction  vet,$t),  one  can  find  the  final  diffraction  point 

(eq,ij>q)  by  following  the  geodesic  path,  step  by  step,  until  the 

geodesic  tangent  coincides  with  (e.,$t).  This  Is  a  rather  tedious  and 
time  consuming  process  If  applied  zfor  each  new  radiation  direction. 
Considering  a  new  radiation  direction  does  not  deviate  greatly  from  the 
previous  direction,  one  should  be  able  to  develop  a  solution  which  uses 
the  properties  of  the  surface  and  the  previous  geodesic  path  to  find 
the  new  diffraction  point.  Such  an  approach  Is  applied  here  to  make 
this  solution  as  efficient  as  possible. 

Since  the  field  decays  exponentially  along  the  ray  path  on  the 
surface.  It  Is  assumed  that  only  one  or  possibly  two  dominant  rays 
exist  In  the  problems  treated. 


B.  SURFACE  GEODESICS 

The  geodesics  on  the  cylinder,  cone  and  spheroid  are  examined  In 
this  section. 

(a)  Cylinder  case: 

The  circular  cylinder  geometry  used  for  this  study  Is  shown  In 
Figure  8(a).  Since  the  cylinder  Is  a  developed  surface,  the  geodesic 
path  Q'Q  Is  a  straight  line  on  the  unfoldedAplanar  surface.  As  shown 
In  Figure  8(b),  the  geodesic  unit  tangent  t  Is  given  by 


A  A  A 

t  *  $COSr  +  zslny  (53) 
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SOURCE 


\lo' 


(b)  UNFOLDED  PLANAR  SURFACE 


Figure  8.  Geodesic  path  on  a  developed  cylinder. 
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where 


a  a  a 

4  a  -sln^x+cos^y. 


(54) 


It  Is  noticed  that  4  and  z  are  the  two  principal  directions  on  the 
cylinder  surface.  For  a  given  geodesic  Q'Q,  one  can  see  that  y  Is  a 
constant  along  the  geodesic  path. 


(b)  Cone  case: 

Consider  a  ray  propagates  along  a  geodesic  path  Q'Q  on  the  cone 
surface  as  shown  In  Figure  9(a)  which  Is  a  straight  line  on  the 
unfolded  planar  surface  as  shown  In  Figure  9(b)  since  the  cone  Is  a 
developed  surface.  Applying  the  same  technique  as  In  Section  II-(E), 
t  and  4  are  the  two  principal  directions  on  the  surface.  The  geodesic 
unit  tangent  t  Is,  then,  represented  by 


A  A  A 

t  =  4cos6+tfis1n$ 

(55) 

where 

A  A  A 

te  *  -xe$1n«+zco$«  , 

(56) 

xe  *  cos4x+s1n4y. 

(57) 

and  6  Is  the  half  cone  angle  as  shown  In  Figure  9(a). 
no  longer  a  constant  along  the  geodesic  path  Q'Q. 

Note  that  a  Is 

(c)  Spheroid  case: 

Recall  that  the  prolate  spheroid  surface  Is  defined  by 
ft(v,4)  *  a  cosvcos4x  +  a  cosvs1n4y  +  c  slnvz 


with 


tanv  ■ 


c  cose 
a  sine 
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Figure  10(a).  Geodesic  path  on  a  sDherold. 


From  Section  II-(E),  It  was  found  that  the  principal  directions  were 

given  by  ^  and  ftv.  As  seen  In  Figure  10(a)  and  (b)  the  unit  vectors 
In  the  principal  directions  can  be  represented  by 


Vi*- 


A  /S 

|  3  -xsln*  +  ycos$ 


tfi  s  Rv/I^vl  =  -x^slna+zcoss 


-a  s1nvxe  +  c  cosvz 

— 5  5  5  9” 
a  sin  v+c  cos  v 


xg  a  cos<j>x+s1n$y 


It  Is  noticed  that  t  Is  merely  4>  .  If  e  denotes  the  angle  between  tf 
and  the  geodesic  unit  tangent  t,  then 


t  *  trcosg+tes1ne  ,  (61) 

which  Is  Identical  to  the  form  used  for  the  cylinder  and  cone  geodesics. 
This  suggests  that  one  might  be  able  to  develop  a  perturbation 
solution  which  gives  a  simplified  form  for  g  on  a  spheroid  using  the 
cylinder  or  cone  expressions  for  &. 


C.  CYLINDER  PERTURBATION 

In  order  to  solve  for  the  geodesics  on  a  prolate  spheroid,  the 
cylinder  perturbation  technique  Is  used  when  the  source  Is  located  at 
e  =90°.  As  shown  In  Figure  11,  the  cylinder  Is  perturbed  by  bending  It 
afound  the  spheroid,  and  the  geodesics  are  to  be  resolved  on  this 
perturbed  cylinder  model.  Note  that  the  source  position  Is  assumed  to 
be  In  the  $$=0  plane.  In  that  different  locations  can  be  dealt  with 

by  simply  rotating  the  pattern  axes. 

Recalling  that  v  Is  a  constant  along  a  given  geodesic  path  Q'Q  on 
the  perturbed  cylinder  In  Figure  11,  one  obtains  a  geodesic  equation 
given  by 


(62) 


tdOy  a  £— 
r 

with  S  aa$.  The  elliptic  cross-section  In  the  x  -z  plane  Is  described 
by  r  e 

x  *  a  cosv' 
e 

z  *  c  slnv' 


with  -x<y<ir.  and  -|<$,  .  At  the  diffraction  point  Q,  the  radiation 

direction  should  coincide  with  the  geodesic  tangent  t  given  In 

Equation  (61).  Thus, 

A  A  A  A 

t  ■  xs1netcos$t+ys1nets1n$t+zcoset 
-  trcosY+tes1nr  (64) 

Equating  the  x-,  y-,  and  z-components  respectively,  one  obtains 
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t  :  -slnacosv- ■■  s^nvC0-s*s^n*  »  s1netcos<|>,. 

12 2  2  2  1  * 
ya  sin  v+c  cos  v 


(65) 


t  :  cos^cosy  -  .  s,„9(s1n4t 

J  fo  o  o  *  t 


r?  ?  ? 

fa  sin  v+c  cosv 


(66) 


t2  : 


C  CQSvSlny 


/ a2s1n2v+c2cos2v 


=  cose. 


(67) 


Now,  t  cos$  -  t  s1n$  yields 

y  ^ 


cosy  =  sine^slnf^-#) 


(68) 


Equation  (67)  Is  divided  by  (68)  to  give 


c  cosvtany  _  c0*  gt 

~2  5  5  2~  s1n (♦*“♦) 

fa  sin  v+c  cos  v 


By  combining  Equation  (63)  with  the  above  expression,  ope  obtains  the 
following: 


f(et,*tv,$)  =  c  cosvs1n($t-$)J0^a2s1n2v'+c2cos2v'dv' 
-a$cotet-/a2s1n^ 


2  2  2 
~  v+C  COS  V 


=  0 


(69) 


Next,  Equations  (67)  and  (68)  are  squared,  which  gives 
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o  o  o  o  ?  2  2  2 

c  wos  vsi n  y  =  cos  et(a  sin  v+  c  cos  v) 

2  2  2 
cos  Y  =  sin  etsin 


From  the  above  two  expressions,  one  obtains 

O2cos2v  =  cos2et(a2sin2v+c2cos2v)  +  c2cos2vsin2etsin2(<t>t-<!s)  . 


Another  function  is,  then,  defined  by 


2  2 
c  cos  v 


n . 


(70) 


Provided  that  one  has  obtained  a  diffraction  point  (v,<(.)  for  a 
radiation  direction  (0t,4>t),  a  numerical  technique  can  now  be  developed 

from  Equations  (69)  and  (70)  to  solve  for  (v+Av,4>+A<i>)  associated  with  a 
new  radiation  direction  ( e t+A0 1 ,4>t+A«t>t) .  Assuming  that  the  ith  set  of 

(et,*t,v,$)  Is  first  known  to  satisfy  f^g^O,  or  at  least 

approximately  so,  the  next  set  (e^+Aej-.^+A^.v+Av.^+Ae-)  Is  obtained  by 

enforcing  f^+^=g^+i=0,  such  that 

fl+r  V\AV%tAVfvAv+V*  =  0  ’  and 
9l+l  =  V\AV\AVVv+V*  =  0  • 

In  matrix  form,  it  is  given  by 
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(71) 


- 

r  1 

f  f 

'  Av 

V  <p 

%  % 

A*  I 

— 

L  j 

K  -  Vet  ■  v*t ! 

i  t  yt  j 

-g,  -  9sMt  -  gtHtJ 


Note  that  the  partial  derivatives  are  given  by  the 

V  n  n  2  2 

f  -  c  r.ir.( 4, r~sin\>J  ya  sin  v'+c  cos  v'dv' 


2  2  2  2 
+  cosv-  jja  sin  v+c  cos  v  ] 


|4.(c2-a2)sin2v.cotet//a2sln2v+c2l 


0 

COS~\> 


f,  =  -C  COZ'-C OS 


-'t 


j|a2sin2v'+c2cos2v'dv' 


I  ?  2  2  2 

-a  /a  si n  v+c  cos  v . cots 


2  2  2  2 
f0  =  a4'csc0t-lia  sin  v+c  cos  v 


v  r  2  o  2  2 

f  =  c  cosvcos  (<}>,.  ••<{>)  J  ya  sin  v'+c  cos  v'dv' 
*t  0 

gv  =  sin2v-[a2cos2et+c2sin2etcos2($t-<t>)] 

2  2  7 

g^  =  -c  cos  vsin  etsin2(<j)t-<t>} 

2  2  2  2  2 

gQ  =  -sin2et[a  sin  v+c  cos  vcos  Ut-4>)] 

2  2  2 

gx  =  c  cos  vsin  o«.sin?U  ,.-4>) 

'♦t  t  t 


Therefore,  one  obtains  (av,A< t>)  for  a  known  ( A6 1 ,A<{>. 


following: 


),  using 
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Equation  (71).  In  order  to  determine  the  initial  diffraction  point 
(v»$)  for  a  radiation  direction  one  can  always  assume  a 

diffraction  point  at  the  source  (v,<(>)  a  (0,0)  with  the  radiation 

direction  (et,±^),  and  gradually  add  increments  (A6t,A4>t)  until  the 

desired  radiation  direction  (et,<|>t)  is  reached.  After  the  initial 

diffraction  point  Is  Identified  by  (v,<|>)}  y>  and  therefore,  the 
geodesic  path  is  determined  by  Equation  (63).  Such  a  numerical 
approach  is  Illustrated  In  Figure  11.  Note  that  one  need  not  trace  out 
the  complete  geodesic  path  from  the  source  location  to  the  diffraction 
point  for  each  new  radiation  direction.  As  shown  in  Figure  11,  the 
diffraction  point  (v+m>,<}>+A$)  for  the  next  radiation  direction  Is 
determined  from  (v,<|>),  using  Equation  (71),  if  ( Ae^. ,A<t>t )  are  small 

which  is  the  case  when  a  complete  radiation  pattern  is  computed. 

After  the  geodesic  path  is  determined,  various  other  parameters 
associated  with  Equations  (15)-(18)  must  be  found.  The  Fock  parameter 
S  of  Equation  (12)  Is  given  by 


Q  ,  kP  1/3 
C  =  /  (-/)  d* 

Q,po  2 


I 


where  the  integral  is  evaluated  along  the  geodesic  path.  Note  that 
i  Is  the  arc  length  along  the  geodesic  and  given  by 


i  - 


slny  CCSy 


Therefore , 


da  = 


s'inY 


/a2sin2v+c2cos2vdv 


or. 


dfc 


J _ 

COSY 


ad$ 


9 


and  the  integration  can  be  rewritten  as 


i 
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(72) 


or 


(73) 


where  Pg  =  l/U^osV^sin^)  *  and  Ki  >Kz  are  9iven  1n  Equations  (45) 


Next,  the  ray  divergence  factor  j defined  as  the  change 

in  the  separation  of  adjacent  surface  rays,  and  It  can  be  found  from  a 
numerical  approach  by  setting  ^constant.  The  ray  geometry  Is  shown  in 
Figure  12. 


Figure  12.  Illustration  of  the  spread 
factor  (/d^/cfp  )  terms. 
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Now,  Equation  (62)  Is  rewritten  as 


cotY 


therefore. 


— 9  ^7  ’i/f  dS  a  -CSC  dy 

uWr7  r 


or 


L  2 

i-*.  ad*  *  -CSC  ydy 

S  J 
5e 


For  a  known  dy=d*0,  usually  set  equal  to  1°,  d*  can,  then,  be 
determined.  Also, 


tan(y+dy)  = 


Sp+/a2si n2v+c2cos2vdv 


Sr+ad* 


The  diffraction  point  (v+dv,*+d*)  of  the  adjacent  ray  Is,  then, 
obtained  from  the  above  equation.  The  radiation  direction  of  the 
adjacent  ray  Is  given  by 


t3  =  t®  COS(y+dy)+ta  Sln(y+dy) 


Recalling  that  the  original  radiation  direction  Is 

A  A  A 

t  «  tpCOSy+tgSlny  , 


(74) 


then  one  obtains 


cos(cty) 


(75) 


I 

In  other  words. 


4, 

W  l7orT(ti‘] 


where  dy  can  be  conveniently  chosen  as  1°. 

This  completes  the  cylinder  perturbation  solution,  which  Is 
applied  If  the  antenna  Is  mounted  on  the  mid-section  of  the  prolate 
spheroid  (e  =*90°,  In  Figure  1  ). 


D.  CONE  PERTURBATION 

When  the  source  Is  not  located  at  the  mid-section  (e  **90°,  in 
Figure  14),  the  cone  perturbation  model  Is  used.  As  shown  In 
Figure  13(a),  the  spheroid  Is  modeled  by  a  perturbed  cone.  The 
associated  unfolded  surface  Is  shown  In  Figure  13(b).  If  y  and  b 
denote  the  angle  between  t  and  t  at  Q'  and  Q  respectively.  It  Is  seen 
In  Section  II I- (B)  that  they  arernot  the  same  as  In  the  cylinder  case. 
In  fact. 


8  *  y-a  (76) 


where 


a  cose 

with  a'  =  acosvs  and  tanv$  =  c  s^ne-  .  With  some  manipulation,  one  can 
show  that  S 

(r-$e)cos(y-a)  =  r  cosy  (77) 


where 


S.  *»  /  /a2s1n2v'+c2cos2v'dv' 
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Geodesic  path  on  the  perturbed  cone. 


rcos  y  * Cr  —  Se )  co*(y  -  a) 


At  the  diffraction  point  (v.$)»  the  radiation  direction  (et»4>t) 
coincides  with  the  geodesic  tangent  t,  and  one  obtains 


t  :  -slnitcosU-g)-  -v-c-0  J s  »  s1netcos<b 

*  r>  2  2 

Ja  sin  v+c  cosv 


t_  :  cos6cos(v-a)-  8  V*y*±W™te&l  =  S1ne*s1n^ 

y  r 2  2  2  2~ 

fj a  sin  v+c  cos  v 


t  .  c  cosvsInlY-gJ.  .  cose 

2  !  2  2  2  2  T 
(a  sin  v+c  cos  v 


Now,  tyCos$-txs1n(j>  yields 


cos(y-o)  =  slno^sl  n(<j>t-4>) 


Using  ( txcos«|»+tysi n<|>)/tz  gives 


-  |tanv  =  tane^cos^^-^)  . 


Three  functions  can,  then,  be  constructed  as  follows  from  Equations 
( 77)- ( 79) : 


f( v,o,y)  *  (r-Se)cos(Y-«)-rcosr  *  0  (80) 

g(e*,4>t,a,Y)  a  cos(Y-a)-s1nets1n(<|.t-  {t  ■*)  *  0  (81) 

h(et,*t,v,a)  *  «  cosets1nv  +  c  s1notcosv  cosUt  -  ~r  a)  (82) 


Using  the  Increment  method  again,  one  uses 
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f1+1  =  f1+fvAv+faAa+fYA^  '  0  » 

91+1  ■  9^+9e  A0t+9^  ^t+gaAa+gyAY  =  0  .and 
t  t 


In  matrix  form,  it  is  given  by 


r f ,  f  fi 

! 

Av 

r~fi  i 

v  a  y 

i 

i 

! 

0  %  *y 

Aa 

-grgetAet‘g$tA*t 

A  ha  0  . 

*Y_ 

-”h1"hetAet"h*tA*t. 

(83) 


where  the  partial  derivatives  are  evaluated  using  the  last  radiation 
direction  results  (et,<|>t,v,<|>)  and  are  given  as  follows: 

l~2  ?  2  2 

fv  *  -Ja  sin  v+c  cos  vCOs(Y-a) 
fa  ■  (r-Se)sin(Y-a) 


fY  =  -(r-Se)sin(Y-a)+rsinY 


ga  =  sin(Y-a)  +  Jr  s1netcosUt  -  Jr  ot) 


=  -s in( Y~a) 

=  -cosetsin(4>t  -  Jr  a) 

=  -s1netcos(4»t  -  jr  a) 

*  a  cose^cosv-c  sine^sinvcosU 


t 


a) 


! 


J 

f 
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h<*  *  c£rSlnetcosvs1n($t  -  jp-a) 

hQ  =  -a  sinets1nv  +  c  cosetcos\>cosU^  -  a) 
t 

h  =  -c  sine  cosvsInU*  -  n-o) 

+t  a 


It  Is  seen  that  one  can  solve  for  Uv.Aa.Ay),  for  a  known  (&et,A4>t) , 

using  Equation  (83).  To  obtain  a  diffraction  point  (v,<j>)  for  a 
radiation  direction  (0*.$*.).  one  can  always  assume  the  first 
diffraction  point  Is  at  tne  source  (v ,<(>)=( v$ ,0)  with  the  radiation 

direction  {£,  ±£),  and  gradually  add  the  Increments  ( Ae^, until  the 

final  radiation  direction  (0t,<l>t)  Is  reached.  Again,  one  need  not 

start  out  from  the  source  everytime,  but  obtains  the  new  solution 
directly  from  Equation  (83),  provided  that  the  new  radiation  direction 
does  not  deviate  greatly  from  the  previous  direction. 

The  geodesic  arc  length  Is  obtained  from  either  of  the  following 
equations: 


fccosy  =  (r-Se)s1na  ,  or 
tsiny  =  r-(r-Se)cosa 


(84) 


Therefore, 


d*  =>  -r-.yO_s.Y  .  da  ,  or 
cos^(y-o) 

_  _  -  __ 
d*  =  jMj.  v<cps_v  dv 
sln(Y-a) 


and  the  Fock  parameter  z  Is  obtained  by  Integrating  along  v  ora,  l.e.. 


K  •  r  cosy/  -L  (-^S.) 


3_ , 


o  pg  L  COS^(y-a') 


da'  ,  or 


(85) 
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e  -  / 


0  pg 


kpG  V3 
(-/)  • 


a2sin2v'+c2cos2\ 


Sln(y-a') 


uv 


(86) 


d^.CQ* ) 

To  evaluate  the  ray  divergence  factor  .  an  adjacent 

surface  ray  is  to  be  determined  as  done  In  the  previous  section. 
Differentiating  Equation  (84)  with  *=constant,  one  obtains 

l~ o  5  2  2~ 

tsinydy  =  sina/a  sin  v+c  cos  vdv  -  (r-Se)cosada, 

I  ~2  2  2  2~ 

tcosydy  =  cosai/a  sin  v+c  cos  vdv  +  (r-Se)sinada,  or 


/  2  2  2  2 
s1no\/a  sin  v+c  cos  v 

-  (r-Se)cosu 

dv 

tsiny 

1  ?  2  2  2 
cosoja  sin  v+c  cos  v 

(r-Se)sina 

da 

fccosy 

For  a  given  dy=d<ji  ,  which  is  usually  set  equal  to  1°,  one  can  find  the 
diffraction  point°(v+dv,<|i+d$)  of  the  adjacent  ray  using  the  above 
equation.  The  radiation  direction  of  the  adjacent  ray  is  then  given  by 

ta  =  tacos(8+dB)+tasin(&+d8) 


where  b  =  y-«  and  dB  =  dy-da.  With  the  original  radiation  direction 
given  by 


t  =  trcos8+tes1nB  ,  then 
>«(  «‘>  r  dr’ 


This  completes  the  cone  perturbation  model  used  to  compute  the 
geodesic  ray  paths  on  the  prolate  spheroid  provided  the  source  Is  not 
mounted  on  the  mid-section  of  the  prolate  spheroid. 
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CHAPTER  IV 


RESULTS 


The  solutions  presented  in  the  previous  chapters  are  employed 
here  to  compute  the  far  field  radiation  patterns  for  slots  or  monopoles 
mounted  on  a  prolate  spheroid.  Various  conical  patterns  are  calculated 
and  compared  with  the  measurements. 

To  show  the  validity  of  the  cylinder  perturbation  solution,  the 
source  is  placed  at  e  =90°  as  shown  in  Figure  14.  To  examine  different 
conical  pattern  cuts,  a  cartesian  coordinate  system  x'-y'-z'  originally 
defining  the  spheroid  geometry  in  Figure  14  Is  now  rotated  Into  x-y-z 
as  shown  <n  Figure  15.  Note  that  the  new  cartesian  coordinates  are 
found  by  first  rotating  about  the  z'  axis  the  angle  <f>  and  then  about 
the  y-axls  the  angle  e  .  The  pattern  Is,  then,  takefi  In  the  x-y-z 
coordinate  system  with  ce  fixed  and  $  varied. 

Calculated  and  measured  azimuth-plane  (e  e90°,  <t>  *0°,  e=90°) 
patterns  for  a  short  monopole  and  a  circumferential  slot  antenna 
located  at  e  =90°,(j>  =0°  on  a  2xx4x  spheroid  are  shown  In  Figure  16(a) 
and  (b).  Here(e  ,  $  )  denotes  the  source  location  In  the  original 
x'-y'-z'  coordinate  system.  More  results  with  ec=90  ,  $  e0°,  but 
different  values  of  e  are  shown  in  Flqures  16(c)  and  (d).  It  Is  noted 
that,  when  e>90°,  the  patterns  are  actually  taken  entirely  in  the 
shadow  region. 

The  geodesic  paths  and  tangents  associated  with  this  source  posi¬ 
tion  are  shown  in  Figure  17(a)(b)  and  compared  with  the  more  elaborate 
work  done  by  Burnside  [2].  Note  that  the  predicted  geodesic  paths 
coincide  with  each  other  In  the  near  region  close  to  the  source. 
Actually, the  Fock  functions  associated  with  the  solutions  drop  more 
than  20  db,  as  the  Fock  parameter  (?)  reaches  2.5  In  the  deep  shadow 
region  [6].  This  clearly  shows  the  slonlflcant  portion  of  the 
surface  as  discussed  In  Chapter  I. 

Several  measurements  were  taken  for  a  short  monopole  placed  at 
o  *30°,  ♦.=0°  to  compare  with  the  cone  perturbation  solution.  Calcula¬ 
ted  and  measured  azimuth-plane  (e  =90°,  $  =0°,  0*90°)  patterns  are 
shown  In  Figure  18.  The  results  are  seencto  agree  very  well  with  each 
other.  The  calculated  geodesic  paths  and  tangents  for  this  case  are 
shown  In  Figure  19(a) (b) . 

Next,  a  roll  plane  pattern  (e  =0°,  ♦  =0° 9  ee90°)  was  examined. 

The  calculated  and  measured  resultscfor  E„  and  EJ  are  shown  In 

e  4* 
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SOURCE 


Figure  20,  where  only  one  ray  Is  employed  In  the  calculation.  However, 
the  one-dominant-ray  assumption  Is  not  adequate  In  this  case,  as  there 
actually  exists  four  significant  6T0  terms  on  a  spheroid  [2]  as  shown 
in  Figure  21.  The  single  exponentially  decaying  field  contributed  from 
either  ray  1  or  ray  2  is  shown  In  Figure  22.  If  one  Includes  ray  1 
and  2  In  the  solution,  one  obtains  the  characteristic  Interference 
pattern  shown  in  Figure  23.  One  would  expect  that  the  further  Inclusion 
of  the  third  and  even  the  fourth  ray  shown  In  Figure  21  might  make  the 
results  agree  better.  However,  ray  3  and  4  fall  Into  the  so-called 
caustic  region  where  virtually  an  infinite  set  of  rays  have  significant 
effects  on  the  pattern,  and  the  basic  GTD  theory  falls.  In  Figure 
19(b),  one  can  see  that  rays  just  simply  collapse  in  these  caustic 
regions. 


The  Ee  pattern  in  the  elevation  plane  (et=90° 
essentially  negligible;  whereas,  the  dominant  E<|>  1 

Here  the  caustic  effect,  which  is  associated  with 


$r*90u,  e=90°)  Is 
n  in  Figure  24. 

is  suppressed 

artificially  so  that  one  does  not  observe  the  sharp  spikes  normally 
attributed  to  a  caustics.  The  discontinuities  In  the  pattern  result 
from  the  fact  that  various  ray  trajectories  are  Included  in  different 
sectors  of  the  pattern.  As  shown  in  Figure  25,  to  remove  the  ray 
tracing  problems  close  to  the  caustic  regions,  a  cone  Is  used  as  a 
boundary  that  defines  whether  two  rays  (ray  1  and  2)  are  to  be  included 
in  the  solution.  It  is  understood  that  the  results  do  not  match  in 
some  regions,  however.  In  that  the  caustic  effect  still  remains 
unresolved. 


To  avoid  the  caustic  effect,  the  pattern  is  taken  again  with  the 
spheroid  rotated  by  45°(i.e.,  e  =90°,  4>c=45°,  e=90°).  The  results 
shown  in  Figure  26  are  again  se6n  to  agfee  very  well  even  with  only  one 
dominant  ray  considered. 

One  has  to  be  reminded  that  the  results  shown  above  are  all  done 
for  a  very  stringent  case,  i.e.  the  spheroid  Is  small  electrically  and 
the  spheroid  is  approaching  a  sphere.  To  model  an  aircraft  or  missile 
fuselage  In  a  physical  sense,  the  electric  dimensions  of  a  spheroid 
would  be  much  larger  than  the  ones  chosen,  i.e.,  2xx4\.  It  Is  noticed 
that  the  source  location  0s=3O°  was  also  chosen  too  close  to  the  end  of 
the  spheroid.  If  more  realistic  aircraft  or  missile  fuselage 
dimensions  were  employed,  one  would  expect  that  the  results  would  be 
improved  a  great  deal.  On  the  other  hand,  it  Is  useful  to  make  such 
comparisons  between  measured  and  calculated  results.  Although 
discrepancies  do  occur  in  some  cases,  they  only  happen  when  the  caustic 
effects  come  Into  play.  However,  the  caustic  effects  associated  with 
this  type  of  caustic  still  remains  an  unsolved  problem. 


More  calculated  examples  are  shown  in  Figure  27-35.  The  algorithm 
boundaries  shown  in  Fi^iffre  25  are  used  in  determining  whether  one  or 
two  rays  are  used  In  the  solution.  Note  that  Bj?  Is  defined  automati¬ 
cally  by  the  numerical  solution.  This  Is  done  'Lby  determining  the 
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caustic  angle  In  the  elevation  pattern  (e  )  and  adding  a  few  additional 
degrees  to  that  value,  l.e.  e^f+AG  whefe  2°<ab<10°.  One  would 
expect  to  observe  slight  discontinuities  somewhere,  because  various 
numbers  of  rays  are  Included  in  different  regions.  Note  that  these 
results  represent  much  more  realistic  fuselage  shapes  and  dimensions. 


<£  (DEGREES  ) 


Figure  16(a).  Azimuth  plane  patterns  (o  =90°,  4>  =0°. 

v  v 

0=90°)  for  a  short  monopole  mounted  at 
e  =90°  on  a  2xx4x  spheroid. 


-  MEASURED 

- CALCULATED 
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£ ( DEGREES ) 


Figure  16(b).  Azimuth  plane  patterns  (ec=90r,  4>c*»  , 

e»90°)  for  a  circumferential  slot  (A=0.291X 
B=0.624X)  mounted  at  fls*90°  on  a  2xx4x 
spheroid. 


<f>  (  DEGREES  ) 


Figure  16(c).  Conical  patterns  (oc=90°,  <f>c=0°,  o=100 

for  a  circumferential  slot  (A=0.291x, 
B=0.624x)  mounted  at  e  =90°  on  a 
2Xx4x  spheroid. 


<t>  (  DEGREES  ) 


Figure  16(d).  Conical  patterns  (E^)  for  ec=90°,  ♦c“0° 

for  a  circumferential  slot  (A*0, 291 x, 
B*0.624x)  mounted  at  e$*90°  on  a 

2\x4x  spheroid. 


53 


(DEGREES) 


PERTURBATION 
REF.  (2) 
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?q  (DEGREES) 

Figure  17(a).  Geodesic  path  defined  by  the  surface  parameters 
(6q,4q),  for  a  source  mounted  at  es"90°  on  a 

2Xx4x  spheroid.  Note  that  y  Is  the  angle  between 
the  geodesic  tangent  €  and  one  principal  direction 
at  the  source  location. 
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Figure  17(b).  Geodesic  tangent  defined  by  the  radial  vector 
direction  for  a  source  mounted  at 

es*90°  on  a  2xx4x  spheroid.  Note  that  y  Is 

the  angle  between  the  geodesic  tanqent  t  and 
one  principal  direction  tr  at  the  source  location. 
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(DEGREES) 

Figure  19(a).  Geodesic  path  defined  by  the  surface  oarameters 
(6q,4>q)  for  a  source  mounted  at  0S*3O°  on  a 

2\xA\  spheroid.  Note  that  y  is  the  angle 
between  the  geodes  1 cotangent  t  and  one 
principal  direction  t  at  the  source  location. 
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<f>i  (  DEGREES  ) 


Figure  19(b).  Geodesic  tangent  defined  by  the  radial  vector 
direction  for  a  source  mounted  at 

0c=3O°  on  a  2xx4x  spheroid.  Note  that  y  Is 
the  angle  between  the  geodesic  tangent  t  and 
one  principal  direction  tr  at  the  source 
location. 


DIFFRACTION 

POINTS 


Figure  21 , 


/ 


The  four  dominant  6TD  terms  that  radiate  at 
(9*900,  4-1450). 


Figure  22.  The  single  decaying  ray  for  a  major 
polarization. 
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TANGENT  PLANE 


TO  SOURCE 


Figure  25.  Cone  boundary  used  to  define  terms  to  be 
Included  In  the  shadow  region. 
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Figure  26(a).  Conical  patterns  (E9)  for  @c*90  ,  ♦c*45°,  and  0»9O  for  a  short 

monopole  mounted  at  0S»3O°  on  a  2\x4x  spheroid.  Only  one  dominant 
ray  Is  employed  In  the  calculation. 


Figure  27.  Radiation  patterns  for  a  short  monopole  mounted 
at  e  =30°  on  a  2,vx4x  spheroid. 


Figure  28.  Radiation  patterns  for  a  short  monopole  mounted 
at  e  *600  on  a  2xx4x  spheroid. 


(d)  ec*90°,  *c=90°,  9=90°  (e)  ec=90°,  4>c=0°,  e*90° 

Figure  31.  Radiation  patterns  for  a  short  monopole  mounted 
at  es=60°  on  a  2xxl0x  spheroid. 


(d)  0  =90°,  $  =90°,  0=9OC 


(e)  0C=9OU,  4>c=0  ,  0=90 


Figure  34.  Radiation  patterns  for  a  short  monopole  mounted 
at  es=60°on  a  2xx50x  spheroid. 
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CHAPTER  V 


SUMMARY  AND  CONCLUSIONS 


The  object  of  this  study  has  been  to  develop  an  efficient 
numerical  solution  for  the  high  frequency  radiation  patterns  of  a 
spheroid-mounted  antenna.  The  Geometrical  Theory  of  Diffraction  is 
used  in  this  study  to  calculate  the  radiation  patterns.  The  curved 
surface  diffraction  solutions  are  discussed  in  Chapter  II,  where  the 
creeping  wave  solutions  In  the  shadow  region  are  of  particular  Interest. 

Cylinder  and  cone  perturbation  methods  are  presented  in 
Chapter  III  to  simulate  the  geodesic  paths  on  a  spheroid,  which  In  turn 
can  be  used  to  model  an  aircraft  or  missile  fuselage.  Because  the 
cylinder  and  cone  are  developed  surfaces,  the  geodesic  paths  can  be 
found  on  the  unfolded  planar  surface.  For  a  given  radiation  direction 
in  the  shadow  region,  the  geodesic  path  and  final  diffraction  point  on 
the  spheroid  can,  then,  be  found  via  an  efficient  numerical  approach. 

The  Geometrical  Theory  of  Diffraction  is  combined  with  the 
perturbation  methods  In  Chapter  IV  to  obtain  various  radiation  patterns. 
The  measured  and  calculated  results  are  seen  to  agree  very  well  except 
when  a  caustic  comes  into  play.  However,  this  type  of  caustic  effect  is 
still  an  unsolved  problem  and  is  beyond  the  scope  of  the  present  study. 

These  solutions  are  very  useful  In  predicting  the  high  frequency 
radiation  patterns  for  antennas  mounted  on  a  spheroid.  The  cylinder 
and  cone  perturbation  techniques  are  highly  efficient  in  approximating 
the  geodesic  paths  on  the  spheroid  for  this  problem.  Further,  this 
perturbation  method  can  be  extended  to  more  general  convex  surfaces 
such  as  ellipsoids. 
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